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We develop a general theory of adiabatic output coupling from trapped, weakly-interacting, 
atomic Bose-Einstein Condensates at finite temperatures. For weak coupling, the output rate from 
the condensate and the excited levels in the trap settles in a time proportional to the inverse of the 
spectral width of the coupling to the output modes. We discuss the properties of the output atoms 
in the quasi-steady-state where the population inside the trap is not appreciably depleted. We show 
how the composition of the output beam, containing the condensate and the thermal component, 
may be controlled by changing the frequency of the output coupling lasers. This composition 
determines the first and second order coherence of the output beam. We discuss the changes in 
the composition of the Bose gas left in the trap and show how non-resonant output coupling can 
stimulate either the evaporation of the thermal excitations in the trap or the growth of the non- 
thermal excitations, when pairs of correlated atoms leave the condensate. 



I. INTRODUCTION 



Trapped atomic Bose-Einstein Condensates (BEC) are now routinely produced in various laboratories around the 
world, and it is important to understand the factors that influence the coherence of atoms transferred from them. This 
is an essential issue for the atom laser research, which has the long-term goal of producing continuous, directional, and 
coherent beams of atoms. A matter-wave pulse was first produced by using a radio frequency (RF) electromagnetic 
pulse to transfer atoms out of a trap, where they were allowed to fall freely under gravity [0-|3j. More recently, 
a stimulated Raman process induced a transition to an untrappcd magnetic state in an experiment at NIST; net 
momentum kick provided by the process resulted in a highly directional beam [Q. On the other hand, a long beam 
of atoms falling under gravity was produced in Munich using an RF-field- induced transition It is noted that 
although the more general features of the output in these experiments are fairly well-understood, detailed properties 
of the atoms in the output beam, and the evolution of the component that remains inside the trap have not so far 
been investigated. 

Previous theoretical treatments of the output couplers for condensates have been either limited to a single-mode 
non-interacting trapped condensate ||-|^ or to mean field treatment for the condensate [|9[-p^ which assume that 
the output beam is extracted out of a condensate at zero temperature and that they can be described by a single 
complex function of space and time. However, real condensates appear at finite temperatures, and as a result, thermal 
excitations play a major role. 

In a previous paper |15|, we outlined a theory of weak output coupling from a partially condensed, trapped Bose 
gas at finite temperatures. By applying the self-consistent Hartree-Fock-Bogoliubov (HFB) theory for Bose gases 
at finite temperatures, we identified three output components. The first is the output of pure condensates, which 
we have called "coherent output." The second is the fraction emerging from the thermal excitations in the trap, 
which are coupled out of the trap by the process of "stimulated quantum evaporation." This is equivalent to the 
quantum evaporation of Helium atoms from the surface of superfluid "^He, where phonon excitations travel up to 
the surface of the superfluid and then spontaneously emerge from the surface as evaporated atoms p^ . In our case 
such an evaporation from the trap is stimulated by an electromagnetic field. The last output component comes from 
the process of "pair breaking," which involves simultaneous creation of an output coupled atom and an elementary 
excitation (quasi-particle) within the trap. For suitable choices of the coupling parameters each of the three processes 
can become the dominant process. We have shown that output coupling can serve not only as a useful way to extract 
an atomic beam out of a trap, but also as a probe to the delicate features of the quantum state of the Bose gas, 
including the pair correlations inside the condensate. 

In this paper, we present an extensive analysis of the spectrum of the output atoms, and address issues that were 
not included in our shorter work. The first is the conditions for the output coupling to give a steady flow of atoms. 
We discuss the behaviour of the output rate and the atomic density in the short and the long time regimes, and also 
discuss the conditions for achieving a steady output beam. Second, the application of output coupling must cause 
changes in the state of the trapped Bose gas, such as changes of the number of excitations relative to the number of 
condensate atoms in the trap. We present a thorough discussion of these changes. The state of the Bose gas in the 
trap is usually described by the Bogoliubov formalism, which assumes an indefinite number of atoms in the system 
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and therefore is not number-conserving. In this paper we discuss a number-conserving description of the system, 
which is especially useful when we consider the process of pair-breaking. 

The results of this paper are directly applicable for any output coupling scheme which involves a single trapped 
state and one output state. We demonstrate here general fundamental issues by considering a one-dimensional Bose 
gas in a harmonic potential, which is coupled into a free output level in the absence of gravity. 

The structure of this paper is as follows: We begin by deriving the equations of motion for the evolution of dynamical 
variables inside and outside the trap in Section O. We present in Section [II a quasi-steady-state formalism, which 
enables us to obtain the properties of the output atoms. We demonstrate the results by a numerical example. We 
outline in Section IV the solution to the equations of motion in the trap that were derived in Section ||, by introducing 
a number-conserving, time dependent Hartree-Fock-Bogoliubov (HFB) formulation in an adiabatic approximation. 
Applying this, we obtain expressions for the internal modes of the system in two different regimes, from which the 
time dependent quasiparticle excitations can be calculated. Finally, discussions and summary are given in Section V. 



II. THE TWO-STATE OUTPUT COUPLING MODEL 

In this section we present our model for describing the output coupling of a trapped Bose gas into free output 
modes. We derive the equations of motion for the atomic field operators in the trapped and the untrapped states, 
and give a general form of their solutions. 

A. Description of the model 

Our model assumes that atoms are initially in an atomic magnetic level \t) ("the trapped state") confined by a 
potential and in thermal equilibrium. A coupling interaction is then switched on, inducing transitions to a different 
magnetic level |/) (the "free" or "untrapped" state). We stress that these are labels denoting internal atomic levels, 
not the centre of mass states, so that for short enough times an atom in an |/) state may still be present within the 
trap. We use the atomic field operator '0t(r) to describe the amplitude for the annihilation of a trapped atom at point 
r, and the operator V'/(r) to describe the corresponding amplitude for a free, untrapped atom. The Hamiltonian of 
the system takes the form 



H — Ti-f^ + + Hcouple, 



(1) 



where 'hif'^ and 'H'^j' describe the dynamics of the trapped and untrapped atoms respectively while Ticoupio 
the coupling between the two states. 

The dynamics inside the trap are given by the many-body Hamiltonian: 



/(O) 



describes 



1 



2y72 



2m 



Vt{v) 



Mr) 



(2) 



where m is the mass of a single atom, Vt (r) is the potential responsible for the confinement of the atoms in the trap 
and Utt is the inter-particle potential between the trapped atoms. 

With the output atoms, significant effect of their (elastic) collisions with the trapped atoms must be taken into 
account. In addition, wc consider a small rate of output from the trap, so that the output atoms are dilute; this 
enables us to neglect the interactions between the free atoms themselves. Since the Bose gases are typically so dilute 
that their mean-free-path for the inelastic collisions is much larger than the dimensions of the atomic cloud in the 
trap, one may neglect also any inelastic collisions of the output atoms with the trapped atoms ijl^. The Hamiltonian 
for the output atoms is then given by 
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Vf{v) 



V'/(r) 



d^v'ij]{v)i;l{v')Utf{v - r')V/(r')V't(r), 
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where V/(r) is the potential that influences the propagation of the output free atoms and Utf is the colhsional 
interaction between the trapped and free atoms; this is in general different from the interaction Uu- We use the usual 
^-function form for the inter-particle potentials 

Uttir - r') = UoS{r - r') (4) 
Utf{r-r') = UiS{r-r'), (5) 

where Uq = 'inh^ att / m and Ui ~ AixYi' atj jm are proportional to the s-wave scattering lengths au and atj for 
trapped-trapped and trapped-free collisions, respectively. We assume a repulsive interaction between the atoms, i.e. 
C/o,i > 0. 

For the Hamiltonian Ticoupio, we consider coupling by an electromagnetic (EM) field which induces transitions 
between the states \t) and |/). In the rotating wave approximation, the EM coupling mechanism may be described 
by the following Hamiltonian, which can quite clearly be generalised to describe any kind of linear coupling such as 
weak tunnelling: 

Wcoupic = ^ y d^rA(r, t)i,) {v)i,t (r) + h.c. (6) 

Here A(r,t) denotes the amplitude of coupling between the trapped and the untrapped magnetic states. The form of 
A(r,t) depends on the type of coupling used: Typical mechanisms are direct (one-photon) radio-frequency transition 
and indirect (two-photon) stimulated Raman transition. For any EM induced processes, the coupling can be written 
as 

A(r,i) = A(r,i)e*('^™-"--^™*\ (7) 

where A is slowly varying in space and time. A can be either time-independent, to describe a continuous electromagnetic 
wave, or pulsed. ?ikom and ?iAeni measure the net momentum and energy transfer from the EM field to an output 
atom. In an RF coupling scheme, A(r, i) is the Rabi frequency il(r, i) = {■p)£{v,t) /Ti (or {jl)B{r,t)/h) corresponding 
to the flipping of the atomic electric (or magnetic) dipole (p) (or (/t)) in the electric (or magnetic) field f (r, t) {B{r, t)); 
Acm is the detuning of the EM field frequency from the transition frequency and kcm is, in general, negligible compared 
to the initial momentum distribution of the atoms. In the stimulated Raman coupling, two laser beams are used to 
induce a transition from |<) to |/) through an intermediate level \i), and 

A(r,0 = "-*'(-'f^'(-'^ (8) 

where Vlu and Qfi are the Rabi frequencies corresponding to the intermediate transitions and A^ is their detuning 
from resonance with the two beams. Acm and kcm are the differences between the frequencies and momenta associated 
with the two laser beams: 

Aem=WlL-W2L T —, (9) 

n 

kem kiL - k2L, (10) 

where eI'^^ — E^p is the energy splitting between the atomic levels \t) and |/) in the centre of the trap. A more 
detailed derivation of Eq. (^ for the Raman process is provided in Appendix An energy level diagram depicting 
the output coupling through the stimulated Raman process is given in Fig. |l]. 



B. Equations of motion 



The coupled equations of motion for the trapped and free field operators are obtained by computing their commu- 
tation relations with the Hamiltonian (Q). We thus find 

-V't(r,i) = -'-Ct^t{v.t) - '-Uoi>\{v,t)Mv,t)Mv.t) 

-zA*(r,t)V3/(r,t), (11) 
^^/(r,i) = -J/:/V^/(r,t)-*A(r,i)^t(r,i), (12) 
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where 



C 



f 



(13) 
(14) 



In Eq. (0) wc have used a mean-field approximation for the colhsional effect of the trapped atoms on the untrapped 
ones. This approximation neglects inelastic scattering processes with the trapped atoms |17 and other possible effects 
of entanglement of the output atoms with the atoms in the trap. The approximation is justified under the assumption 
of long atomic mean-free-path mentioned above. 



C. General solutions 



1. Output atoms 
The formal solution of Eq. (|l2|) for ^/ in terms of is 

i^f{v,t) = i;f(r,t)-ij^ dt' J £v'Kf{vy,t-t')\{v',t')M'r',t'), (15) 

where ip'^p satisfies the time-dependent Schrodinger equation for the free evolution of tjjf in the absence of output 
coupling 

and the "free" propagator Kf{r, r' ,t — t') satisfies the partial differential equation 

^^Kf{r, r', t-t') = -^CfKfir, r', i - t') + <5(r - r')S{t - t'). (17) 

The second term in Eq. ( p^ describes the transition of the atoms from the trapped level \t) into the free level |/) 
with amplitude A(r, t) and their subsequent propagation as free atoms. 

We note that it is useful to expand the field operator ipf in terms of the normal modes </5k(r) of the untrapped level 

Vi/(r,i)=^^k(r)6k(i), (18) 

k 

where &k satisfy the usual bosonic commutation relations [^k,^k/] = ^k.k'- k denotes the momentum state of the 
free atoms with energy = hLO\^, and (/Sk are the time-independent solutions of the single-particle problem in the 
non-trapping effective potential V/(r) + Ui{ipl{r)ijjt{r)) created by the mean field effect of the collisions between the 
trapped and the free atoms. In principle, the solutions ip]^ and the energies may change with time due to the 
change in the density of the trapped atoms and the subsequent change in the effective repulsive potential near the 
trap. In what follows we will neglect this time-dependence under the assumption of weak output coupling and slow 
changes in the density of the trapped atoms. 

In the absence of gravitational or other forces, at positions far away from the trap, k may be taken to be the wave 
number of a plane wave (/3k ~ e*'' '' with Wk — hk^ /2m. In the presence of gravity, the structure of the output modes 
should be defined appropriately, and the modes k may be given asymptotically by the solutions of the Schrodinger 
equation in a homogeneous field. In Eq. (^8|) we have used a sum over discrete output states. It is noted that the 
actual structure of the Hilbert space for the output modes depends on the potential V/(r); if this potential vanishes 
far away from the centre of the trap, then the sum should be replaced by an integral / d^k and the operators 6k 
should be defined such that [Sk, bl^,] = 5(k — k'). 

In terms of the basis functions iy9k(r), the free field operator is given by 

V-f (r,t)=^¥.k(r)6k(0)e— * (19) 

k 

and the propagator of the free atoms may be written as 
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Kf{r,r',t-t') 



k 



(20) 



It is useful to describe the evolution of the untrapped atoms in terms of the annihilation operators 6k in a specific 
mode k. The solution for this operator is obtained by multiplying Eq. (p^) by ip^ and integrating over all space: 



6k(i) = 6k(0)e 



dt' / d3r(^^(r)A(r,t')e-*"''(*-*')^t(r,i')- 



(21) 



Solutions for the output field '0t(r, i) in Eqs. (15) and ( pl| ) require an explicit expression for the trapped field 
operator ipt{r,t'). The simplest approximation is to take the first order solution in the coupling amplitude A(r,t). 

This corresponds to a very weak couphng and ipt (r , i) ~ il^^^ (r, i) , where ip^^^ (r , t) is the field operator of the trapped 
atoms without output coupling. The fundamental properties of the output under such approximation will be the main 
subject of Sec. III. 



2. Trapped atoms 



By substituting the solution ( |15D for ipf back into Eq. ( |ll| ) we obtain the following equation for the trapped field 
operator V't(r, t) 



-7^t(r,i) = ^^-CtM^,t) - ^Uoi^Ur,t)Mr,t)My^,t) 



dt' I d\G(v,v' ,t,t')M^' ,t') - i\*{v,t)i^f {Y,t), 



where 



G(r, r', t, t') = A*(r, t)Kf{v, r' , t - t')X{r', t') 



(22) 



(23) 



The solution of the integro-differential equation, Eq. (^2|) will be the main subject of Sec. |rV|. In principle, two 
different situations may be expected from such an integro-differential equation. In the case where this equation 
describes a coupling to the output levels with a narrow available bandwidth compared to the coupling strength A(r, t), 
we anticipate Rabi oscillations of the atomic population between the trapped and untrapped levels. Physically, this 
means that the output atoms stay near the trap for a long enough time to perform these oscillations. On the other 
hand, when the bandwidth of the output modes is large compared to the coupling strength, an exponential decay of 
the population in the trap is expected. Physically, this behaviour is expected when the output atoms are fast enough 
to escape from the trap before the interaction couples them back into the trapped level. Even if the coupling is very 
weak, an oscillatory kind of behaviour is expected for short times compared to the inverse of the bandwidth of the 
relevant output modes, before the output rate settles on a constant rate with fixed energy. For this last case, Eq. (12? 
may be viewed as a Langevin equation for an interaction of a confined system with an infinite heat bath [ p2[ . 

The trapped field operators if^ti^^t) can, in principle, be expanded similarly to the field operator ipf{r,t) as 



^t(i",i) = ^<?!'n(r)a„(i). 



(24) 



where 0„ (r) are the normal modes of the trap given by the solutions of the single-particle problem in the potential 
well Vt{r). These eigenmodes, however, do not form a good basis, since the interaction between the atoms results in 
a strong mixing between the levels. In the following sections, therefore, we will instead use the basis of condensate 
and its quasi-particle excitations, which are obtained from the Hartree-Fock-Bogoliubov theory of an interacting Bose 
gas. 



III. OUTPUT PROPERTIES IN THE QUASI-STEADY STATE 



In this section we present the properties of the output atoms under the assumption that the output coupling is 
very weak. In this case, the output beam of atoms can serve as a probe to the structure of the Bose gas in the trap 
under steady-state conditions. The solution for the output is given by Eq. (p^), or equivalently Eq. pl|), where we 
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substitute V't(r, t) « 'ip'f\r, t). We start this section by giving a brief description of the formahsm that allows the use 
of the steady-state solution ij)f\r,t) for the atoms in the trap. We then present basic properties of the output, in 
particular the spectrum and the density of the output coupled atoms. Finally the first and second order correlations 
of the output atoms are presented. 



A. The trapped atoms in the quasi steady-state 

We briefly review in this subsection the theory of the trapped Bose gas in steady-state conditions, in a way that 



will enable us in Sec. IV to extend the theory to the time-dependent case where the number of atoms in the trap 
does change adiabatically during output coupling. Since the situations discussed in this paper involve transitions 
of atoms into untrapped propagating states where counting the number of output atoms could be one of the main 
measurements, we choose to use a number-conserving theory in the spirit of the theories put forward recently pO| . In 
addition, the theory enables an extension of the finite temperature HFB-Popov method into time dependent cases. 

For describing a partially condensed system of atoms at a finite temperature, we split the field operator for the 
atoms in the trap into a part which is proportional to the condensate wave function and a part which represents 
excitations orthogonal to this state: 

V-tlr, t) = e-''*(*) [^o(r)ao W + 5i^{v, t)\ . (25) 
Here $(t) is a global phase given by 

$(0= f Kt')dt', (26) 
Jo 

where /i, the chemical potential of the system for the given global variables, is constant under steady-state conditions. 
The operator oq is a bosonic annihilation operator satisfying [ao,aJ] = 1 and describes the annihilation of one atom 
in the condensate state ipo{r). The number of condensate atoms is represented by the operator Nq = ajao. The 
non-condensate part, 6ip, is assumed to be orthogonal to the condensate in the sense J d^ri(;Q{r, t)5il}{Y, t) = 0. In the 
number conserving formalism 5il) is approximated by the following Bogoliubov form: 

<5V'(r,i) « ao^^K-(r)a,(t) +t;;(r)a](i)], (27) 

where dj, d] are the bosonic operators satisfying [di, dj] = 5ij in the space of states with non-zero condensate number, 
and they describe the annihilation or creation of excitations (quasi-particles), or, equivalently, transitions from an 
excited state j into the condensate and vice versa. This implies that the operators a.j , d| do not commute with the 
condensate operator clq. The wave functions Wj(r) and Uj(r) are the corresponding amplitudes associated with the 
annihilation of a real particle at position r, an action which involves both annihilation or creation of excitations on 
top of the condensate. The time-dependence of the functions '4)Q,Uj,Vj is induced only by the change in the global 
variables Vt{v,t), Nt(t), Etrap{t) and they are assumed to be time-independent under the steady-state conditions. 

The condensate wave function V'o(r) in Eq. ( p5| ) is defined as the solution of the generalised steady-state Gross- 
Pitaevskii equation 

{Ct~^i + [/o[^o|^o(r)P + 2n(r)]}7Ao(r) = (28) 

where Ct is given in Eq. (^3|), while the adiabatic mean number Nq of the condensate atoms and the density n{v) of 
the non-condensate atoms are calculated self-consistently by requiring 

No + j d^rn{r) = Nt. (29) 

The functions Uj{r),Vj{r) satisfy the steady-state equations 

Ct -^i + 2Uo[No\Mr)? + n{r)] C/o^o(^o(r))2 W w,(r) \ / w,(r) 

-C/oA^o(^o(r))' -Ct+^i- 2C/o[iVo|Vo(r)|2 + n(r)] J { v,ir) J ' { v,{r) 

+U0N0 J d^r\Mr)nro{r)u,{r) + Mr)v,{r)] ) , (30) 
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where Ej = huj are the jth quasi-particle excitation energy with respect to the condensate ground state energy, and 
the second term on the right hand side ensures the orthogonahty of the non-condensate functions with the condensate 

pl[ . We note that the vectors ^ ^ satisfy an equation similar to Eq. ( |30| ) with Ej —Ej; Eqs. ( psf ) and ( |30| ) 

must be solved self-consistently for any given global conditions. The mean number of atoms in an excited state in 
equilibrium is assumed to be given by the Bose-Einstein distribution 

" eMfi^i/T]-l' ^^^^ 

In this section we assume a weak output process so that the total number of atoms in the trap does not change 
significantly during the application of the output coupling. Under these conditions the time-dependence of the 
operators ao and dj is assumed to be simply 

ao{t) w ao(0) (32) 
aj{t) ^ dj{0)e-''^'K (33) 

In our numerical demonstration throughout this paper we take a one dimensional Bose gas of Nt — 2000 atoms 
in a harmonic trap with frequency lu. The critical temperature for condensation in this case is Tc ~ 300?iw/fc. We 
have used a self-consistent HFB-Popov method to find t he wave functions and energies of the condensate and the 
excitations. Throughout this paper we take Uo = IQhuj^fm] moj. We present calculations for two temperatures: for 
T — lOhuj/k we obtain the chemical potential /z ~ 2.5huj and the non-condensate fraction ~ 2%. At T = 150'huj/k we 
obtain /i « 2.3?ioj/k and the non-condensate fraction is ~ 44%. We take the interaction strength between the trapped 
and untrapped atoms to be Ui — Uq- Length will be presented in units of 2y^h/muj ("harmonic-oscillator units"). 



B. Basic properties of the output 

In order to obtain the properties of the output we expand the output field operator in terms of the free modes 
(/?k in the quasi-steady-state [Eq. (p^]. We assume A(r, i) — A(r)e~'^='"*. Using the form ( p5| ) and ( |27| ) of tpt and the 
assumptions (^2|) and (js^), we obtain from Eq. (|2^) the following equation for the annihilation operators of the free 
output modes 



6k(i) « e-*""* {6k(0) - z {AkoCko(t)ao(0) 



-ao-L ^[Akj+i^k,+(t)«,(0) + Ak,_i?kj- (0^1(0)] \ \ , (34) 

V ^0 7 



where 



Ako = J d3r(^* (r)A(r)Vo(r) (35) 
Ak,+ = J d\^l{r)Xir)u,{r) (36) 
Ak,- = / dW^ir)X{r)v*ir) (37) 



are the matrix elements of A(r) between the wave functions of the collective excitations and the output states. The 
time- and energy- dependence is determined by the functions 

e"*('^out~"k)* — 1 

Dkr, it) = I , (38) 

^out ~ 

where 

=M + Ae,„+^^, (39) 
with 77 = Ojj+jj — , and Ej+ = Ej and Ej^ = —Ej. 
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With the above definitions, the field operator tpf{r,t) of the free atoms can be written as 
V'/(r,t)«V;f (r,i)-**°(r,t)ao(0) 



-ao(0)- 



J2 (r, t)a, (0) + ^- (r, t)a] (0) 



(40) 



where 



*/(^'^) =II'^k(r)I?kr,(i)Akr,e 



(41) 



for each 7] ~ 0,j+, j—. 

This resuh enables us to calculate various properties of the output from the trap, assuming Bose-Einstein statistics 
for the initial quasiparticle populations inside the trap. For the initial state we assume that all the cross correlation 
fimctions between different operators aQ,aj,a'j vanish, and the only non-zero contributions are the populations 



{a]aj) = 



(42) 
(43) 
(44) 



Any measurable quantities related to the output atoms may be expressed in terms of the correlation functions of 
the field operator ?/;/(r,t) at different times and space points. In particular, the density of output atoms at a given 
point r and time t is given by the equal-time, equal-position correlation function 



nout{r,t) = {ipj{r,t)-ipf{r,t)). 



(45) 



By using Eq. ( [40| ) and the assumptions (42)-([44|) we observe that ?T.o„t(r,t) can be written as a sum over discrete 
contributions from the levels in the trap 



nout{.r,t) = 7Vo|*/(r,Or +1] (r,t)r + (", + 1)1^-}" (r,i)p 

j 



(46) 



where the first term is the condensate output, the second term is the contribution from the stimulated quantum 
evaporation of the thermal excitations, and the third term is the contribution from the pair-breaking process, as 
discussed in Ref. p5| . 

Following similar steps, the number of output atoms in mode k of the free atomic level, rik = {blJjk), is given by a 
sum of discrete contributions from the different levels of the trapped gas 



(47) 



where each term rj — 0,j+,j— in Eq. 



has the form 

Kit) - \Xkr,\^\D^M' 



The time- and energy- dependence of each of the is given by the function 



\Dkvit)\ 



sin[(u;k -t^out)V2] 



(wk - t^o«t)/2 



(48) 



(49) 



This function has a spectral width in Wk which decreases with time as ~ n/t. This spectral width represents the 
energy uncertainty dictated by the finite duration of the output coupling process. The time evolution of the output 
atoms is therefore governed by the spectral dependence of the matrix elements Akr;. 

In order to analyse the evolution of the output rate and the output atom density, we define two frequency scales 
with regard to the matrix elements Ak^ for each 77 — 0,j+,j — : (1) Aw,;, the frequency bandwidth within which the 
matrix elements Akj;, defined in Eqs. (|35|)-(|3^), are significant and (2) Sujri, the "width" of Ak,j in the Wk = hk^/2m 
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space in the vicinity of uj^ = '■o^out- The weak couphng assumption is justified if the strength of the couphng, which 
may be represented by the parameter A = ^ J d^v\\{v)\^ is much smaher than Aw,, of each trap state, namely, 

A < Atj„. (50) 

If this condition is not satisfied, then we expect Rabi oscillations between the trapped atomic level \t) and the output 
level I/} |2^. In the case of weak coupling, we may identify three temporal regimes: 

1. Very short times, t <C Aw"^. Then the function Dkri{t) in Eq. ( ^8|) becomes independent of k, and if uj^ut 
lies within the bandwidth Acj,,, I?kr((i) ~ <• In this case, the completeness of the set of functions 'Pk(r) implies 

^^^{r,t)^t^oMr)ijl'{r)t, (51) 

where V'? = V'o,''/'*^ = ^^'^ i^t^ — The initial shape of the output wave functions before it had time to 
propagate is therefore the overlap of the electromagnetic field amplitude and the corresponding trapped wave 
function. The density noMt(r) in this case is then similar in shape to the density of atoms in the trap and the 
total number of output atoms increases quadratically in time. This result may also be used to calculate the 
output beam immediately after the application of a strong coupling pulse, before the output beam starts to 
propagate or Rabi oscillations occur. 

2. Intermediate times, Atj,y^ < t < Sw^^ ■ In this case the rate of output from each trap state r] may show 
oscillations, which follow from interference between output from different momentum states. 

3. Long times, t 3> 5uj^^. The output from the internal state 77 is then mainly generated in a narrow range of 
energies around fiuj^y^^ and the rate of output dn\^/dt into these specific modes settles on a constant value, which 
is determined by the absolute value of the matrix element Ak,) at Wk — ^^ut- It is then given by the Fermi 
golden rule 

^=2^5]<|Ak„P<5(a;k-^o\*), (52) 

and the output rate obtained when one scans the frequency detuning Ae^ of the coupling fields measures the 
magnitude of the matrix elements |Ak))P as a function of ciJk. 

The asymptotic behaviour of |£'kjj(OP this limit is 

|i?k,(OI' - 2^'5(wk - u^lut)i + 27 ^ir^, (53) 

(Wk -t^o„i) 

where the principal part in the second term is defined as P/ [lo — lo'^Y = d/duj'^[P/ {uj — lo'^)]. The first term represents 
a linearly increasing mono-energetic contribution of the level 77 to the output with a constant rate given in Eq. ([S^). 
The second term involving the principal part represents a time-independent non-rcsonant part that has two physical 
implications. First, it represents those frequency components that are different from the central resonance frequency 
which arise due to the sudden switching on of the output coupling field at t = 0. Second, this term results in the 
fields ^'^(t) containing a non- propagating (bound) part 

(*:;(r,t)),oW = -*P^ ^^i^^, (54) 

which stays mainly near the trap. This term appears as a part of the dressed ground-state of the coupled system, 
which is a mixture of the trapped and the untrapped atomic levels. It therefore represents a virtual transition to the 
output level while the atoms remain bound to the trap. It is detectable if the atomic detecting system is sensitive 
enough to identify small number of atoms in a different Zeeman level near the main atomic cloud, which would 
contain atoms in the Zeeman level \t). Although this last contribution is in general small compared to the resonant 
contributions, it may be significant when considering the condensate output (77 = 0), which is multiplied by the large 
number A'q, and therefore may be dominant near the trap relative to the contributions of the stimulated quantum 
evaporation (for Ag^ < 0) or the pair-breaking (for Aem > 0). The condensate contribution can be estimated by 
{n^f)bound ~ 2A*o(A(0)/(jJo,jj)^, where A(0) is the Rabi frequency associated with the coupling field at the centre of the 
trap. 

The spectral widths Acj^ and 5uJr, defined above may drastically vary with the structure of the Hilbcrt space of the 
output modes, which is determined by the potential V/(r), and also with the spatial shape of the coupling A(r, t). It 
is worth mentioning the following limiting cases: 
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1. In the absence of gravity the wave functions ip^ are roughly given by the plane waves e^^ '^. Then the matrix 
elements Ako that couple the condensate to the free modes is roughly the Fourier transform of the condensate 
wave function V'o(r)- If no momentum kick is provided, then their width in momentum space is given in terms 
of the spatial width tq of the condensate by Sko ^ I/tq. The corresponding spectral width, Aujq, is then 

AujQ ^ Sljq ^ h/2mrQ < ui, (55) 

which implies that the time it will take to achieve a constant rate of output from the condensate is larger than 
the period of the trap. If the condensate is broadened by a strong coUisional repulsion then this time may be 
much greater than this period, which is typically of the order of 10 ms. 

2. In the case where a momentum kick ke,„ is provided, the spectral width for the condensate output becomes 

dujo h\'kem\6kQ/2m ?i|kem|/2TOro. (56) 

This makes the time for achieving steady output shorter by a factor (|kem|ro)^^ compared to the previous case. 
If |kem| corresponds to an optical wavelength then this factor may be of the order of 10. 

3. In the presence of gravity the spectral width of Akr; is determined mainly by the gradient of the gravitational 
potential over the spatial extent of the corresponding wave function ip^. A typical value of this gradient for the 
condensate wave function in the experiment of Ref. [gl is about Sujq 27r x 10 kHz. In this case the time needed 
to achieve a steady output is much shorter, in the order of ^ 0.1 ms. 

The rate of transfer of atoms into the output level as a function of Ae,„ in our one-dimensional example is plotted 
in Figure |^. This rate is a sum of contributions from the condensate and the excited states in the trap 

dt dt ^ 

j 

where 

d^vnl^,{,v)^Y.< (58) 

k 

The rate of output from the condensate component, dn'j/dt (solid line), that from stimulated quantum evaporation, 
dn''i^ /dt (dashed line), and finally from pair-breaking, dn''^^ /dt (dash-dotted line) are shown for temperatures 
T = IQTiLo/k O.OSTc, bold line) and T = IbOTiu/k (~ 0.5Tc, thin line), in the case where no momentum is transferred 
from the EM field (kem — 0). The threshold below which the condensate output is not produced is at Agm = — 
which is slightly different for the two temperatures. To prevent unphysical effects that follow from the divergence of 
the density of states at small momenta in one-dimensional systems, we have assumed that the density of momentum 
states per energy is constant, p(wk) = 1- As to be anticipated, the composition of the output beam changes as a 
function of Aem : for negative values of Aem the dominant contribution is from the stimulated quantum evaporation 
of initially excited levels in the trap; for positive Aem it is the contribution of pair-breaking that dominates the 
output. The output contains overwhelmingly condensate components at central values of Aem- Comparison of the 
results for T = lOtko/k and T — IbOtku/k shows that output rate from the pair-breaking is pronounced mainly at low 
temperatures. 

Figure ^ is a one-dimensional demonstration of the output density for coupling frequencies in (a) the stimulated 
quantum evaporation regime {Aem — ~5uj), (b) the coherent output regime (Aem = 0), and (c) the pair-breaking 
regime {Aem = 8a;). At very short times the output density from each level has the same shape as the density of 
the given level in the trap, as follows from Eq. (^T|). After a short time, the output atoms emerge mainly in two 
momentum states (/3k corresponding to the right- and left-propagating waves with energy hcu^ut ■ Since the magnitude 
of the matrix elements Akjj for these two modes are equal, the output beam corresponding to a given component rj 
forms a standing wave and consequently the density n^^t {x) becomes oscillatory. This aspect is demonstrated below, 
when we discuss the coherence of the output. These standing waves are not expected in cases where the inversion 
symmetry is broken, such as the case with \iem ^ 0. In cases (a) and (c), where Aem is very positive or negative, one 
can see that the output density from the condensate has a steady component that remains near the trap. This part 
corresponds to the appearance of the bound states discussed above in conjunction with Eq. (p^. 



dt 



dt 



(57) 
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C. Coherence of the output 



The concept of the n-th order coherence in a quantum system was originaUy developed in the optical context to 
quantify the correlations in the field . The first order coherence measures the fringe contrast in a typical Young's 
double slit experiment, while the second order coherence gives indications of counting statistics in, say, Hanbury- 
Brown-Twiss experiments. A theory of the coherence of matter-waves was presented only recently p4| and the case 
of a trapped Bose gas has also been discussed [|5|. It follows that, for matter waves, the theory of coherence which 
is applicable to real experiments is much more complicated than that for the optical case. However, any measures 
of coherence must involve correlation functions between the matter-wave field operators. For simplicity, we use here 
definitions of matter wave coherence functions that are equivalent to the optical definitions by replacing the usual 
electromagnetic field operators by the matter-wave field operators. 



1. First-order coherence 



The first order coherence function gj.- {r,r' ,t,t') for the output atoms is defined as 

g) 'ir,r ,t,t) = , ^ / (59) 



where g^^-* = 1 implies full coherence and 5'^^ = implies total incoherence. The first-order coherence for a random 
or thermal mixture of many modes typically takes the maximal value for r = r' (i.e. (7(^-'(r,r) = 1) and falls down to 
zero for large |r — r'| or \t — t'\. Highly monochromatic beams, however, are characterised by the fact that g^'^'' — 1 
even for large |r — r'| or \t — t'\, implying high fringe visibility even if widely separated parts of the beam interfere. 

An atomic beam weakly coupled out from a finite temperature Bose-gas is, in general, a mixture of quasi- 
monochromatic beams originating from the condensate and the internal excitations in the trap. The nature of this 
mixture depends on the frequency, shape and momentum transfer from the electromagnetic field, and correspondingly 
the coherence properties are significantly affected. Following the quasi-steady-state solution for '0/(r) in Eq. ( ^0| ) we 
find for our first order coherence, 

9f\r^r'^t,t')^ ' E N4^Jir,t)r^Jir'^n (60) 

\/nout[r)nout[r ) ,,^oj+j- 

The coherence is maximal if only one of the terms from the sum over 77 is dominant. Fig. ^ shows the first order 
coherence function g^j^\xi, X2,t) of the output atoms in our one-dimensional demonstration as a function of X2 for 
fixed xi ^ and time t = lOO/cj. When Aem ~ and the temperature is low (T — lOhto/k, Fig. ^) the coherence 
fmiction is unity except for points where the condensate density vanishes (see Fig. ||b). At T = IbOTiu/k (Fig. ^) the 
thermal component is larger and it is more dominant near the points where the density of the condensate component 
is low. These features are unique to configurations where the output has a form of standing matter- waves. When 
^em = ^Soj (Fig. I^c) the thermal components are dominant (see Fig. ^) and the coherence drops much lower than 
unity. When Aem = (fig. ^) only few thermal output components exist from the pair-breaking process, and 
consequently one obtains a comparatively high coherence function. 



2. Second-order coherence 



Of particular interest are the intensity correlations of the fields which are important, for example, in experiments 
involving non-resonant light scattering from an atomic gas pq]. These intensity correlations are expressed in terms 
of the second order correlation function, (ri, r2, ti, ^2), which is defined as 

(2), ^ ^ , (V'/(i-i,ii)V'}(r2,i2)V'/(ri,ii)?/'/(r2,t2)) 

Of (ri, ti; r2, t2 ) — — rr . (ol) 

^ (V^}(ri,ti)V'/(ri,ti))(V'}(r2,t2)^/(r2,t2)) 

The function measures the joint probability of detecting two atoms at two space-time points. If the detection prob- 
ability of an atom is independent of the detection probability of another atom then g^^^^ = 1 and the probability 
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distribution is Poissonian. This is the case for a coherent state of the matter field. However, for a thermal state 
the correlation function at the same space-time point is (ri = r2,ti = ^2) ~ 2. This implies that the atoms are 
"bunched," i.e. there is a larger probability to detect two atoms together. 

The second-order correlation function at equal position and time points ri = r2 and <i ~ was previously calculated 
for a trapped Bose gas ||2^. Here we shall follow the same treatment for calculating the second-order coherence of 
the output beam. We decompose the field operator i/)/(r, i) into a part proportional to the condensate and a part 
proportional to the excited states in the trap, and apply Wick's theorem to the expectation value of the product of 
four non-condensate operators: 

(V'L(r)V'L(r)^„c(r)^„c(r)) - 2n^{v) + TO*(r)m(r), (62) 

where we have defined n(r) = (?/'Jj^(r)?/'„c(r)} and TO(r) = (7/'„c(r)^nc(r)). One then obtains for the second order 
coherence of the output atoms: 

gf{v,t) = 1 + {2Re [nO„,(r)n„„,(r) -t- [(vl/0„,(r))*]2m„„,(r)] -HnL(r) + \mout{r)\'} , (63) 

where no„t(r) = EjKJt(r) + niZt{r)] and mout{r) = J2j «'}~(2nj + 1). 

We note that although in Ref. [p7| the terms proportional to rh had negligible contribution, here they may play an 
important role even at zero temperature in situations where the tuning of the coupling EM-field frequency yields an 
output beam that emerges mainly from the non-condensate parts of the trapped gas. 

The equal-time, single-point intensity correlation of the output atoms after a time t — lOO/oj in few typical cases is 
shown in Fig. |^. If the output condensate is dominant (Fig. |^a) the function g^'^\x) is equal to unity except at discrete 
points where the output condensate wave function vanishes. At a higher temperature (Fig. ||b) the thermal output 
components tend to raise g^'^\x) near the points where the coherent part is small. In the case where the thermal 
component is dominant (Aem = — 5w, Fig. ||c) ^'^'(a;) assumes the value of 2. In the case where the pair- breaking is 
dominant (Fig. |^d) the intensity correlations tend to assume values greater than 2. This can be interpreted as an atom 
bunching effect caused by the combination of the process of pair-breaking with the stimulated quantum evaporation 
of the thermal states. 



IV. EVOLUTION OF THE TRAPPED GAS 



The last section was devoted to a discussion of the properties of the output in the quasi-steady-state approximation, 
where the Bose gas in the trap is assumed to remain unchanged by the output coupling process. In this section we 
describe the internal dynamics of the trapped atomic Bose gas during the output coupling process. During this 
process, the trapped atomic population of the condensate state and each of the excited states change in a different 
way and the system is driven out of equilibrium. In the typical case where the duration of the coupling process is 
short compared to the duration of relaxation processes at very low temperatures ||30| , pT|| , the dynamics is represented 
by approximate solutions of Eq. (^2|). In the case of weak coupling, the solutions are best represented in terms of the 
adiabatic basis of the system, which are the steady-state HFB-Popov solutions for a given total number of particles 
and given total energy of the system. It can serve as a good basis as long as the changes in the conditions in the trap 
are slow enough compared with the trap frequency. 

We begin by first introducing a two-component vector formalism that is convenient for dealing with the many 
modes of excitations. We then obtain linear equations of motion for the creation and annihilation operators of the 
condensate and excitations in the adiabatic basis. In the adiabatic conditions these equations may be simplified and 
solved analytically. A perturbation solution is then presented, which is suitable for describing the short time evolution. 
Finally, we find that the number-conserving formalism fails to describe the evolution of the condensate number in the 
pair-breaking regime. This problem is discussed and cured in the end of this section. 



A. Vector formalism for the trapped atomic gas 

The dynamics of the excited states in the trap is usually described by a set of two coupled equations of the form 



Eq. which was discussed in Sec. HI, or its time-dependent version This form, as well as the fact that 

the expansion for the field operator in Eq. ( p5| ) involves the quasiparticle creation and annihilation operators aj , , 
motivates the introduction of the two-component vector formalism as follows. 
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First, we define the normalised condensate operators 

Co = ao— 4 = — F^«o> (64) 

which are well-defined in the space spanned by states with non-zero condensate number; within this space, they satisfy 
cqCq = cjco = 1- We then define the two-component column vector operator 
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ir^t)=( ff;'^':'}) (65) 



which describes transitions from the condensate state to itself and to and from the excited states. 

The expansion of V't(r, t) in Eq. ( p5| ) and is equivalent to the expansion of (r, t) in terms of the two-component 
wave function vectors as 

°° ft 
it{r,t)= ^ C,(r,0«r,(t)e-'Jo'^*^''(*'), (66) 

ri— — oc 

where the index 77 takes integer values from —00 to 00. The index rj = corresponds to the condensate; negative 
indices stand for solutions of Eq. ( ^0| ) with negative energies E^j = —Ej, and operators such that a_j ^ '^]- 
addition, we note that the condensate operator 



do = clcLQ = \/ No (67) 
is Hermitian. The two-component vectors ^,,(r, i) are defined as 

co-fsVc..r!0,^-..r!i ). (68) 



The time dependence of the vectors £,rj{r, t) and the energies Ejj{t) in Eq. ( |66|) is governed by the time-dependence of 
the global variables of the system, while the time dependence of the coefficients d^, rj G (—00, 00), represents changes 
in the populations of the condensate and the excited states. 

The usual orthogonality and normalisation conditions for the eigenfunction uj and vj are written in the vectorised 
notation as 



dh^lir)Ur) = J dh^l{r)a3Ur) - 0, (69) 
J d''r^l{v)asUr) - sign{£;j<5^„ (70) 

for any 77, ^ 0. Here 

e] ^ (u* V*); ; ^ (v, u,) (71) 
are the two component row vectors and (73 denotes a 2 x 2 matrix, 



B. Equations of motion for the operators 



We now derive the equations of motion for the operators corresponding to transitions from the condensate to 
the adiabatic eigenmodes of the system and vice versa. We first multiply Eq. ( p^ ) by cje**. The resulting equation, 
together with its Hermitian conjugate, form a set of equations which can be expressed in the following vector form 
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= / dt' / d3r'e»-3[*(t)-*(t')]G(r,r',t,i')6(i-',i') 



where G(r, r',i,t') and A(r) are the matrices 



z~XHr)as^f\r) (73) 



cl{t')G*{r,r',t,t')co{t) 



G(r,r',M')=( ^oK^'^y^,^.^,^J^o,^> , ), (74) 



and in a similar manner to (r, i) of Eq. 



^j°''(r, i) describes the free evolution of the output field operator ipfir, t), as given in Eqs. (|l|), @. The term (^j'"^ 

is the operator describing the free evolution of inside the trap in the absence of the output coupling but with a 
given adiabatic change in the global variables. Here we use the same approximations as in Eqs. (|3^), (|3^), which is 
equivalent to 

{Lr\v,t) = -^Y,E,i,a,e-'^o^^^''^'r (77) 



The time derivative of on the left hand side of Eq. (73) may then be written as 



The first term corresponds to the time dependence due to the change in the global variables, the second term is due 

to the change in the populations of the condensate and excited states, while the third term cancels with (Cf )*-"■' on 
the right-hand side of Eq. ([73|). We multiply Eq. (^, in turn, by ^^0-3 for every i] ^ and by i^o for ?y = 0, and 
integrate over r. This multiplication should be understood as an inner product between row and column vectors. By 
applying the orthogonality and normalisation relations in Eq. (|6^) and Eq. ( |70| ) we obtain the required equation of 
motion 

a,, = -^AV(t)d,(i)-^ J dt'G,j,{t,t')a,{t')-t J rf3rF„(r,t)C';.°)(r,t). (79) 



M,At) = e'/o f d3,^t (r)^,,^,(r) (80) 



Here 



is a matrix with zero diagonal, which describes mixing between the adiabatic levels that is induced by the change in 
the global variables. This term in Eq. (|7^) may be neglected in the adiabatic limit where the change in the global 
variables is very slow. Its effect in slightly non-adiabatic conditions will be discussed elsewhere [32[. 

The second and third terms in Eq. ( [79|) describe changes in the trap which are directly induced by the output 
coupling. Defining 

%{t) = / dt'[^iit') + a3E^it% (81) 

and 
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7? > 

77 = , 
r)<0 



(82) 



one may write 

Gr,At,t') = Jd'rJ d3r'eJ(r,t)a,G(r,r',<,i')e'"^[*''(*)-*-^*')le.(r',i'), (83) 

and 

i^,(r)=^t(r)^^^3e-3*„X*(r) (84) 

which describes the effect of the zero-field fluctuations. An exact analytical solution of Eq. (|7^) is, in general, not 
possible. However, in the following we present two methods of approximate solutions to this equation: an adiabatic 
approximation, which is suitable for describing the evolution at long enough times, and a perturbative expansion, 
which is suitable for short times. 



C. Solution via adiabatic approximation 



First we consider the adiabatic and quasi-continuous case where the functions f,, change very slowly with time 



and the coupling amplitude is given by A(r, t) = A(r)e' 



In this case we let M^i^ « 0, Second, Eq. (|79|) is 



further simplified by finding an approximate expression for the integral involving Grii,{t,t'). We make a Markovian 
approximation, which transforms the integro-differential equation into an ordinary differential equation, which 
can then be solved analytically. Following the definition of G{r, r', t, t') in terms of the free output modes denoted by 
k [Eq. ( pO| ) and (p3|)], the functions Griu{t,t') may be written as 



Grju(t, t') 



i[t:Jk-Ae,„]fT3(t-t')gi[*^(t)-*.(t')])^^^(i')^ 



(85) 



where 



(86) 



with the matrix element Ak,, as defined in Eqs. (^-(^. The time dependence of the matrix elements Ak», is induced 
only by the change in the global variables, which is assumed to be slow. The sum over k in Eq. ( ^ ) may then be 
regarded as a Fourier transform of the products A]{^^CT,,Aki/ over Wk at "fixed" point in time t = t ~ t' . The width 
Ar^y of GrfuiT) as a function of r is then roughly given by the inverse of the spectral width Aa;,,^ of the product 
of the matrix elements Akr; and Akjy, which is, in turn, given by the smallest of the spectral widths Atj,, and Awj^ 
of the corresponding matrix elements. In the same conditions that allow the weak coupling approximations done in 
when G,j^ATjj[y <C 1 and t ^ At,,^, we may take a^(f') « ctjj{t) in Eq. ( [79| ) and cJ(t)co(i') = 1, and 
-oo, namely 



Eq. (|52|), i.e., 

extend the integration over t' to 



it), 



(87) 



where 



Jo 1, 



1 



Wk - M - f^aE^ - Aem - ic^ae 



Ak.(i)- 



The complex fraction should be understood as 



P 

±TrS{x) — i — , 

X 
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where Pjx means the principal part of l/x when integrating over x. 

Further simpUfication is achieved when we notice that if the terms F^j^ are much smaller than the energy splittings 
— Ey between the excitation levels in the trap, then the cross-terms with rj ^ i> oscillate as fast as ~ gi{E^-E^)t 
and their contribution averages to zero. We then obtain a system of separate uncoupled equations for each operator 
djj, which is given, for non-negative ^7 = j > 0, by 



for which the solution is 

cej{t) = exp[- 

For j = 0, 



(89) 



(90) 



(91) 



and for j 7^ 0, 



where 



-3± 



(92) 



(93) 



The imaginary part of Tjj represents energy shifts induced by the output coupling, while its real part 7^ = ReFjj is 
given by 



7j 



± = ReFji = ± ^ |Akj± p(5(t^k - ^JiTE.J- Ae„) 



(94) 



representing decay (7^+ > 0) or growth (7^- < 0) of the population of excited level j. 

We now proceed to calculate the number of condensate and quasi-particle excitations inside the trap under the 
adiabatic approximation. The evolution of the condensate number is given straight-forwardly by 



iVo(i) = (dg) - 7Vo(0)e-27o*. 



(95) 



However, for calculating nj(t) = (aj(i)Q;j(t)) we must also consider the free term in Eq. (90), whose contribution is 
proportional to the correlations of the free field operators (r, t) 



(^f (r,i)(^f )t(r',0)=E^k(r)e^ 



(96) 



In the case of very weak coupling, where Tjj may be assumed to be time-independent, the contribution of this last 
term in Eq. ( |9C| ) to nj{t) is 



1, 



C^2+72 ' 



(97) 



where 



ImTjj. The spectral width of the integrand is Auj ^ n/t for \jjt\ <^ 1 and Aij 



7j 



for 



|7ji| 3> 1. Under the conditions that led to Eq. (|89|), one may take w in the (5-function and consequently identify 
Sk |Akj-P<5(a;k — ^iZt) ~ 7j-- The integration over lu may be then performed to give the final result 
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1 _ e-^-Ti* 

Ujit) = exp[-27ji]nj(0) - 27j_ 

2jj 



(98) 



This equation is the solution of the differential equation 

dji, ■ 

^ = -2jj+n,it) ~ 27,_[n,(t) + 1]. (99) 

Here, the first term on the right-hand-side is responsible for an exponential decrease in the number of excitations due 
to stimulated quantum evaporation, while the second term is responsible for an exponential increase in the number of 
excitations due to the process of pair breaking, which may start even when the excited states are initially unpopulated. 
This increase in the number of excitations must, quite clearly, lead to the increase in the number of atoms in the 
excited states, together with an increase in the number of output atoms. There is, however, no process that may 
balance this growth in the total number of atoms, and this implies that the growth must be compensated by a decrease 
in the number of condensate atoms. This is, in fact, not evident from the above equations and the problem is discussed 
at the end of this section. 



D. Solution via perturbation theory 



A full solution of the linear integro-differential equations Eq. ( [79| ) may be sought by perturbative iterations, taking 
the magnitude of the coupling strength A as a perturbative small parameter. Here we present the second-order 
perturbative solutions, which are valid at short times when the population in different excitation levels are not 
changed significantly from their initial value. In this case we may also assume that the wave functions and energies 
of the condensate and excitations are not changed significantly from their initial values (i.e. M^ii, « 0). 

If we take the zeroth order solution to Eq. ( [79| ) to be given by Eqs. ( ^2[ ) and (33), then the second order solution 
is given by 



a^(0) 



-t dt' J d^rF,ir,t)if\t). 
Under the above assumption, we may perform the integration to obtain 



I k 

-*E^lw^k,(t)(|^|;). 



a.(0) 



Here 



where the functions Dk?; are defined in Eq. (^8|) and 



D 



k?7 

D 



k, — 



= 



l-H^3K~Ae„-M]--E^}t-e-'t'3['^k-Ae^->.]-fi,}t 



By using the identity 

[see Eq. (^9|)] we obtain the following expression for the number of condensate atoms in the trap 



(100) 



No{t) = {alit)) = No{0) 



l-El^kon^ko(t)P 



(101) 



(102) 



(103) 



(104) 



(105) 
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and for the population of the excited levels we obtain 



(106) 



Comparison of Eqs. ( |105| ), (106) with the equivalent expressions for the number of output atoms in Eqs. (pj|), ( |4g 
shows that exactly one condensate particle is taken out of the trap per each output atom generated by the coherent 
output process, while one excitation (quasi-particle) is taken from the trap per each output atom generated by the 
stimulated quantum evaporation, and one excitation (quasi-particle) is created per each atom that leaves the trap 
through the pair-b reaking process. 

From Eq. (101) it is straightforward to compute the correlations {ali{t)a,y{t)) between the condensate and the 
excited levels and between the different levels in the trap. However, it may be shown that only diagonal terms rj = v 
are growing in magnitude with time, while the off-diagonal correlations remain small even after a long time and 
represent the effects of mixing between different levels induced by the coupling interaction. 



E. Number of particles and energy 

The above treatment of the evolution of the system of a trapped Bose gas has used a formalism which conserves 



the total number of particles in the system. However, Eqs. (|95|), (105) show that the change in the number of 
condensate atoms in the system is independent of the changes in the number of quasi-particles in the trap. This leads 
to an apparent violation of number conservation; this violation is most pronounced in the process of pair-breaking in 
which output atoms are created together with quasi-particles in the trap. The problem arises because we ignored the 
off-diagonal part in the Hamiltonian which is also responsible for the changes in the number of condensate atoms, i.e. 

Hoff-diag = UoJ (fr (^o*(i-))' aS4V'nc(r)V'nc(r) + H.c. (107) 

This part of the Hamiltonian, which is responsible for the generation of quantum entanglement between the condensate 
and the excited states, is washed-out in any mean-field treatment such as the HFB-Popov treatment used here. In the 
mean- field theory the time-evolution of the condensate operator ao in the steady-state is simply given by Eq. (|3^), 
and this leads to the apparent violation of number conservation when the number of quasi-particles in the system is 
changing. A rigorous theory which corrects this fault is beyond the scope of this paper. Such a theory is in principle 
straightforward, but technically a little complex: we have to incorporate the anomalous average (V'nc(r)V'nc(r)) into 
the calculation of the condensate wave function and show how this anomalous average acquires an imaginary part in 
the presence of output coupling of excited states. This, from another viewpoint, represents the change in the effective 
T-matrix for the interaction potentials in the presence of decay. Here, we will incorporate number-conservation by 
requiring that the number of condensate atoms iVo(i) is to be determined from the conservation of the total number 
of particles. If the evolution is adiabatic then some time after the switching-on of the coupling interaction the mixing 
between different quasi-particle levels may be neglected. The total number of atoms in the trap is then given by 

N,{t)^N,{t)+Y,L,{t) I dMK(r)p + |t;,.(r)P]+ / d'r\v,ir)A . (108) 
j 

On the other hand, we must require 

Nt{t) ^ NtiO) - Noutit). (109) 



If we compare Eqs. ([47|), ( |4§| ) with Eqs. ( |105| ), ( p_06D we see that in the process of stimulated quantum evaporation 
(r] = j+) the number of quasi-particles in the trap decreases in the same rate as the number of output atoms increases, 
while in the pair-breaking process {rj = j—) the number of quasi-particles in the trap increases in the same rate as 
the number of output atoms increases. In other words, in the stimulated quantum evaporation process one thermal 
quasi-particle is transformed into a real output atom, while in the pai r-br eaking process one quasi-particle is generated 



per each output atom that leaves the trap. From inspection of Eq. (lOS), this implies that for each atom that leaves 
the trap in a stimulated quantum evaporation(SQE) process, the number of particles associated with the quasi-particle 
j in the trap decreases as 
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SNf'^'' = -J dh[\u,{r)\' + \v,{r)\'] = -l-2 J d'r\v,{r)\\ (110) 

This must be compensated by an increase in the condensate atom number by 

SN^'^^ ^ +2 J d^r\v^{r)\\ (111) 

On tlic other hand, in the pair-breaking (PB) process, the number of particles associated with the quasi-particle j in 
the trap increases by 

SNf^ ^^ + ^J '^^i-|wj(r)|2. (112) 
This must be compensated by a decrease in the condensate atom number by 

SN^B " ^'^'■I^J Wl' - -2 (^1 + J d^r\vj{r)\^^ . (113) 



These considerations lead us to corrections to Eq. (^5|) which contains only the changes in the condensate particles 
originating from direct output from the condensate component of the Bose gas. The rate equation for the condensate 
atoms is now 



dt 



-2-iqNq 



2Y,\j d'v\v,{v)\^\^^+ j dMu,{v)\^^Y (114) 



where dn'^'^'^/dt and dnf^/dt are the first and second terms on the right-hand side of Eq. (p9|). The solution of 
Eq. ( |ll4| ) should now replace the previous solution for No{t) in Eq. (p5|). 

The plots of the time evolution of the trapped condensate and non-condensate populations for few temper ature s 



and coupling parameters are given in Fig. |6|. These plots are solutions of the differential equations (|99| ) and ( 114 ). 
When Aem = (Fig. ||a,b) the condensate part decreases while the thermal part does not change significantly. When 
= — 5w (Fig. I^c) conservation of energy only permits transitions from upper excited states to the output level, 
and the population of the condensate and the lower excited states thus remains unchanged. The upper excited 
states, in this case, are found to depopulate completely as to be expected. When A^m — 8uj (Fig. |6|d) the thermal 
population grows significantly due to transitions of unpaired atoms from the condensate into the excited states. 
However, the energy distribution in the lower excited states is a highly non-equilibrium distribution and dissipation 
and thermalization effects that have not been taken into account in this paper should play a major role. The short 
time limit i.e. < t < lOw"^ behaviour is clearly not accurately described in these plots but it can be calculated by 



using the low-order perturbative expansions of Sec. [VD 



Changes in the total energy in the trap may be caused either by the transfer of atoms out of the trap or by the 
changes in the chemical potential n and energies Ej of the excitations. The second kind of process is beyond the 
scope of this paper, since we have neglected changes in /i and Ej and put M^ii, = in Eq. (^9|). As for the first kind 
of process, an energy quantum of 6E — fi leaves the trap for each condensate atom that leaves the trap (consequently 
the energy in the trap is reduced by /i), the energy changes by SE^^^ = — /i — Ej for each atom that leaves the trap 
by the stimulated quantum evaporation process, and finally SEj'^ = — /i + Ej for each atom that leaves the trap 
through the pair-breaking process. Therefore we have a relatively simple result for the rate of change in energy: 

dEt dNf v-^ „ drii 



V. SUMMARY AND DISCUSSION 



In this paper we have set up a general theory of weak output coupling from a trapped Bose-Einstein gas at finite 
temperatures. The formalism developed here is suitable for the discussion of both Radio-frequency or stimulated 
Raman output couplers. It has enabled us to gain much information on the basic features that we expect in real 
experiments: the time-dependence of the output beam, the effects of excitations in the trapped Bose gas and the 
pairing of particles. Predictions for specific systems can also be based on our theory. 
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For the time-dependence of the output beam, we have shown that the output beam is a mixture of components 
from different origins in the trap. The output condensate (ry = 0) is the coherent part of the beam, while each excited 
level j in the trap contributes two partial waves: one originating from the process of stimulated quantum evaporation 
(77 = J-I-), where a quasi-particle (excitation) in the trap transforms into a real output atom, and the other originating 
from the pair-breaking process (77 = j — ), where two correlated atoms in the trap transform into a quasi-particle in 
the trap and a real output atom. We have shown that a steady monochromatic wave from each component is formed 
after a time comparable to the inverse of the bandwidth of the corresponding matrix element Akij as a function of 
Wk. We have also analyzed the oscillatory behaviour of the output rate at short times and showed the existence of 
non-propagating bound states in the untrapped level that are formed near the trap as a result of the mixing induced 
by the output coupler. 

As for the evolution of the Bose gas in the trap during the process of output coupling, we have shown that for 
the case of weak coupling an adiabatic approximation may be made, which enables calculation of the composition 
of the Bose gas inside the trap in terms of the adiabatic basis of condensate and excitations. We have shown that 
exponential decay of the excitations is expected when the stimulated quantum evaporation process is dominant, while 
an exponential growth of the number of excitations is expected when the pair-breaking process is dominant. We 
have shown that the number of trapped condensate atoms increases in each event of stimulated quantum evaporation, 
while it decreases by more than 2 atoms per each event of pair-breaking. However, we stress that a more elaborate 
number-conserving theory of time-dependent evolution of the Bose gas in an open system than that considered here 
is needed. 

The coherence of the output beam was shown to depend on parameters under experimental control such as the 
detuning of the laser. We note that the coherence of the output atoms also tells us about the coherence properties of 
atoms inside the trap; the coherence of the trapped Bose gas is expected to be altered as a direct consequence of the 
output coupling. In simple terms, when the output atoms are mainly those of condensates we expect the coherence 
of the internal atoms to drop, if only because the amount of coherent condensate fraction decreases. The coherence 
of trapped atoms, although interesting theoretically, is not experimentally verifiable. 

Apart from designing an atomic laser with well-controlled beam properties, we saw that the measured output 
properties may be an excellent tool in investigating the nature of trapped Bose gases at finite temperatures. The 
properties of the output beam may be a probe to the temperature of the trapped Bose gas as well as the internal 
structure of the ground state and the excitations. The present treatment may be extended to cope with other possible 
configurations that are likely to appear in the future such as a trap with multi-component condensates and Bose gases 
with negative scattering lengths. 

Finally we note that the pair breaking process, and indeed the output coupling of the condensate in general, provides 
an experimentally feasible method to study quantum entanglement in a macroscopic system. The quantum theory 
of entanglement is currently under intense study owing to its relevance to quantum computation; so far it has rarely 
been studied in the context of BEC. 
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APPENDIX A: EXPRESSION FOR A(R, T) - STIMULATED RAMAN SCHEME 

We derive here the expression for the effective coupling function A(r, t) in Eq. (^) for the stimulated Raman transition 
coupling scheme. A detailed analysis of a Raman coupling process from a condensate in the mean-field approach can 
be found in Ref. 0. 

We consider a single atom which can be found either in the trapped level \t) or in the free level |/}. A pair of laser 
beams with spatial and temporal amplitudes EtL(r,t) and EfL{r,t) are responsible for non-resonant transitions from 
\t) and I/) to a high energy level \i). The Hamiltonian is then given by 

H= J2 [?i^,+fff]U)01 + ^ E[/^r.^^-^(^'^)l*)01 + /^-C-] (AI) 

Here htOj are the internal energies of the levels and 
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where Vj (r) are effective external potentials acting on the different atomic levels, and fijs are the dipole moments for the 
transition |j) — > The amplitudes of the laser fields are assumed to have the form EjL{r,t) = fj(r, t)e*'''j^ '"~"j^*' 
where the envelopes SjL are slowly varying with respect to the exponential term. The wave function describing the 
atom has the form X!j=t / i The equations of motion for the three amplitudes are: 

^h^t = V* - fi:,E:^^Pre-^'-'''-'^^^' (A3) 
tn^Pj = - /i},i?;iV.e-^(— (A4) 

ihi;, = V» - f^uEtL^te-'^^^*-^'^' ~ fifLEfL^fe-'^^f-^^'>\ (A5) 



The solution of the equation for -0^ as a function of the two other amplitudes can be written as 



h^°^(r f)^i f dt' f d\'K,{r,r',t-t') V e-'("^-"')*' 
X (r', i') V'i (!■', t')e'^^^^ -r'-c^,!,*') . (A6) 



Here -i/'l"-' is the solution of the Schrodinger equation ihdipi/dt — Hf'^'ipi with the initial condition -0'^*'^(O) — ipiiO) — 0, 
under the assumption that level \i) is initially unpopulated. U,j{r,t) = fiji£jL{'r,t)/2h are the slowly varying Rabi 
frequencies and Ki is the propagator for the evolution of the level \i), which can be expanded in terms of the energy 

eigenfunctions of in a similar way to the expansion in Eq. (pO| 



The crucial step now is to notice that the main time-dependence in the time-integral in Eq. (AC) comes from the 
terms Q-^i^i-'^i+^iL)t ^ whose frequency of oscillation is assumed to be in the optical range, while the other terms, 
which correspond to atomic centre-of-mass motion are oscillating in frequencies below the radio-frequency range. 
Assuming that the switching time of the coupling is much longer than the short period of oscillation of the fast terms, 
we expect the contribution to the integral in t' to come only from a short time interval around the end-point t. We 
then take t' —t'va the slow terms. As a result, we have i^i(r, r', t — t') ~ Ki{r, r', 0) = 6{r — r') and t') ~ V'i(r, t), 

and obtain 

0,(r,t) = - ^ r!,(r,i)^,(r,i)e*'^^-'-. (A7) 

UJj — LOi + UJl 



When this is substituted in Eqs. (|A3|) , ( |A4| ), and the rapidly oscillating terms are dropped, we obtain 



= -J^Hf^^t + «Att(r, t)4>t{v, t) + iXtfiv, t)^f{r, t) (A8) 
= -J^^/°V/+*A//(r,i)0/(r,i)+zA/t(r,t)V't(r,t), (A9) 



where 



UJy - -I- Wj/L 

The form of A(r.t) in Eq. (|^) is achieved by assuming iot — LOi-V lOtL ujf — uii + ojfL = and then noticing that 
Xtf = Ajj. We have also neglected the diagonal terms Xu,Xff, which are responsible for an additional effective 
potential acting on the levels |t) and |/), under the assumption that they are small compared to the other potentials 
Vt(r) and V/(r) near the trap. This assumption is justified in the adiabatic case discussed in this paper, where the 
coupling is assumed to be weak and slow. 
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FIG. 1. Schematic diagram of energy levels and couplings involved in the stimulated Raman process. A is the detuning, oji, QiL, 
i = 1,2 are the frequencies and the Rabi frequencies of the two lasers. 
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FIG. 2. The rate of output as a function of Aem for atoms emerging from various different processes at temperatures T = lOhui/k 
(bold line) and T = 150?ia)/fc (thin line). Solid line: dn^j/dt, output from the condensate component; Dashed line: y^^. dn?^ /dt, from 

stimulated quantum evaporation; Dash-dotted line: dn'j~ /dt, from pair-breaking. A constant density of states p(u;k) = 1 was used. 
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FIG. 3. A one-dimensional demonstration of the temporal evolution of the coherent component (bold line) and the thermal component 
(thin line) of the output atomic density at T = IbOhui/k 0.5Tc) for different coupling strengths and detunings. (a) Aem = —5a;, 
A = 0.5a;, (b) Aem = 0, A = 0.2u!, and (c) Aem = 8w, A = 2u!. The output density from the condensate has a steady component that 
remains near the trap; this part corresponds to the appearance of the bound states discussed after Eq. fe4h. 
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FIG. 4. The first order coherence gj. {xi,X2,t) of the output atoms as a function of x = X2 for a fixed value oi xi = at time 
t = lOO/tJ, for (a) T = 10hLj/k,Aem = (Dominant coherent output), (b) T = 150hw/k, Aem = 0, (c) T = 150?itj/fc, Aem = -5uj 
(dominant thermal output) and (d) T = I50hu>/k, Aem = 8t<; (dominant pair- breaking). 
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FIG. 5. Equal-time single-point second order coherence function g^^\x) of the output atoms at time t = 100/i.iJ. Figures (a)-(d) 
correspond to the same cases plotted in figure n of the first-order coherence. 
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FIG. 6. Evolution of the condensate (solid thin line) and excited (dashed line) atomic populations in the trap for few temperatures 
and coupling parameters. The bold solid line shows the total number of atoms in the trap as a function of time. In (a) (T = lOhLo/k) 
and (b) (T = 150hui/k) the output is dominantly from the condensate. In (c) (Aem = —5a;) stimulated quantum evaporation from the 
higher excited levels is dominant and the remaining population in the trap is mainly the condensate and the lower excited levels. In (d) 
(Aem = So)) the population of the lower excited levels increases due to pair-breaking, while the condensate depopulates. 
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